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Abstract. We present a newly measured X-ray temperature function of galaxy clusters using a complete flux- 
limited sample of 61 clusters. The sample is constructed with the total survey area of 8.14 steradians and the 
flux limit of 1.99 x 10“^^ ergs s“^ cm“^ in the 0.1-2.4 keV band. X-ray temperatures and fluxes of the sample 
clusters were accurately measured with ASCA and ROSAT data. The derived temperature function covers an 
unprecedentedly wide temperature range of 1.4-11 keV. By fitting these data with theoretically predicted tem¬ 
perature functions given by the Press-Schechter formalism together with a recent formation approximation and 
the CDM power spectrum, we obtained tight and individual constraints on Dm,o and erg. We also employed the 
Formation-Epoch model in which the distribution in the formation epoch of clusters as well as the temperature 
evolution are taken into account, showing significantly different results. Systematics caused by the uncertainty in 
the mass-temperature relation are studied and found to be as large as the statistical errors. 
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1. Introduction 

The mass function of clusters of galaxies (MF), the num¬ 
ber density of the most massive virialized systems, con¬ 
tains information on the structure formation history of 
the universe. A theoretical framework, e.g. the Press- 
Schechter formalism (Press & Schechter 1974) together 
with the Cold Dark Matter model, has been established 
to predict the MF. This allows us to constrain cosmo¬ 
logical parameters using an observationally determined 
MF for the present epoch (as well as its time evolution 
for even tighter constraints). In particular, ag, the ampli¬ 
tude of mass density fluctuations on a scale of Mpe 
where h is the Hubble constant in unit of 100 km/s/Mpc, 
and the mean matter density, are most sensitively 

determined by the cluster abundance measurements (e.g. 
Henry & Arnaud 1991; White et al. 1993| ; Eke et al. 1996| ; 
Kitayama & Suto 1996; Viana & Fiddle 1996; Oukbir & 


Blanchard 1997; Pen 1998; Eke et al. 1998). 


The estimated virial masses for individual clusters depend 
rather strongly o n mod el assumptions, however. As argued 
by Evrard et al. (1997) on the basis of hydrodynamical N- 
body simulations, cluster masses may be presently more 
accurately determined from a temperature measurement 
and a mass-temperature relation determined from detailed 
observations or numerical modeling. Thus alternatively, as 
a well-defined observational quantity, the X-ray temper¬ 
ature function (XTF) has been measured, which can be 
converted to the ME by means of the mass-temperature 
relation. 

The first measurements of the XTE were reported by 
Edge et al. (199C) and Henry & Arnaud (1991), using 
an X-ray flux-limited sample of 45 clusters and 25 clus¬ 
ters, respectively. Recent observational improvement on 
the XTF was made by Markevitch (1998), Henry ( ^000 ), 
Blanchard et al. ( ^000 ), and Pierpaoli et al. (2001), using 
more accurate temperature-measurement results for each 


Observationally the local MF has been derived from 
measuring masses of individual clusters from galaxy veloc¬ 
ity dispersions or other optical properties by Bahcall and 
Cen ( 1995 ), Biviano et al. ( 1993D , and Girardi et al. ( 19981) . 
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cluster with ASCA data (Tanaka et al. 1994). However, 
the narrow temperature ranges (3-10 keV) in which the 
XTF is defined so far do not allow an investigation of a 
more detailed shape of the XTF other than just fitting a 
single power-law. Without better information on the ac¬ 
tual shape of the XTF, no independent constraints on erg 
and Dm 0 can be derived, and a wide range of combinations 
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is still allowed. Recently, a number of new X-ray-cluster 
surveys were performed, which provide a high complete¬ 
ness for the brightest clusters. This motivated us to re¬ 
visit the measurement of the local XTF and to improve 
its accuracy in order to derive narrower constraints on the 
cosmological parameters. 

Reiprich & Bohringer ( ^002|) have compiled a new X- 
ray flux-limited cluster sample (HIFLUGCS) with a flux 
limit of 2 X ergs s“^ cm“^ (0.1-2.4 keV). Compared 

with previous samples with similar flux limits, their cat¬ 
alog covers the largest volume and is the most complete. 
Based on this cluster sample, Reiprich & Bohringer have 
measured total masses for individual clusters and derived 
the X-ray mass function for the first time. In this paper we, 
using this cluster sample, report the construction of a new 
measurement of the XTF. The larger number of clusters in 
the sample together with accurate temperature measure¬ 
ments with ASCA data leads to a significant improvement 
of the XTF, which covers 1.4-11 keV temperature range. 
The X-ray flux-limited sample and the new temperature 
measurements with ASCA data are presented in Sect. 2. 
The derivation of the XTF is described in Sect. 3. In Sect. 
4, the XTF is used to obtain constraints on cosmologi¬ 
cal parameters, and the results are compared with other 
work in Sect. 5. Throughout the paper, the Hubble con¬ 
stant is given as 100 h km/s/Mpc, and log and In denotes 
a decimal and natural logarithm, respectively. 


2. Sample and analysis of ASCA data 


2.1. Master catalog 

For this study we use a sample of 106 g alaxy clusters com¬ 
piled by Reiprich & Bohringer ( 2002|) , which is used to 
construct a flux-limited complete sample (HIFLUGCS) 
with the flux-limit of 2.0 x 10“^^ ergs s“^ cm“^ (0.1- 
2.4 keV). The 106 clusters have been derived from combin¬ 
ing previous X-ray cluster or elliptical galaxy catalo gs tha t 
include Bohringer et al. ( 200C , ^001 ), Bohringer ( 1999 ), 
Retzlaff et al. in prepa ration, Ebeling et al. ( 1996 , 199S ), 
de Gr and! et al. ( 1999 ), Beuing et al. ( 1999 ), Lahav et al. 
( 1989 ), and Edge et al. ( 1990| ). Most of them are based 
on the ROSAT All Sky Survey (Triimper 1993 ; Voges et 
al. 1999 ). From these catalogs, all objects satisfying cer¬ 
tain criteria, mainly that the measured flux given in any 


of these catalogs is above 1.7 x 10 ergs s 


( 0 . 1 - 


2.4 keV), have been sampled. A defined region around 
the center of M87 is not included in the survey area, be¬ 
cause of the large scale emission of the Virgo cluster that 
compromises cluster detection and characterization in the 
region. This excludes the Virgo cluster and M86 from the 
master catalog. The count rate of the PSPC (Pfeffermann 
et al. 1987) in the 0.1-2.4 keV band for each object was 
redetermined from the ROSAT All Sky Survey data or 
PSPC pointing observations whenever available by inte¬ 
grating all X-ray flux within a maximum radius where the 
cluster emi ssion is significantly detected. See Reiprich & 
Bohringer ( 2002| ) for detailed descriptions of the sample 


construction and data analysis. The names and the red- 
shifts compiled from the recent literature of all the 106 
clusters are listed in Table 1. 

As described below, we newly determine the X-ray 
temperature as well as the X-ray flux for each of the 106 
clusters. ASCA data are available for 88 of the 106 clus¬ 
ters, and the X-ray temperatures are measured by analyz¬ 
ing the ASCA data (Sect. 2.2). Using the ASCA results, 
the correlation between the X-ray luminosity and the tem¬ 
perature, L — T relation, is established (Sect. 2.3). Among 
the other clusters for which no ASCA data exist (18 clus¬ 
ters) , temperatures of 4 clusters have been measured with 
Einstein, EXOSAT, or XMM-Newton in previous works 
and were used in our analysis. For the remaining 14 clus¬ 
ters, the temperatures are estimated from the observed 
ROSAT PSPC count rate by means of the L — T rela¬ 
tion. By setting a flux-limit of 1.99 x 10“^^ ergs s“^ cm“^ 
(0.1-2.4 keV), a complete flux-limited sample comprising 
63 clusters is then constructed (Sect. 2.4). 


2.2. Temperature Measurements 


For our study we need to determine an average tempera¬ 
ture for each cluster that closely reflects the mass of the 
cluster, i.e. a virial temperature. In many clusters, the X- 
ray emitting hot gas is not isothermal but a cooler gas 
component is often observed in the central region. The 
cool component, which is due to either a cooling flow (e.g. 
Fabian 1994) or the ISM of cD galaxies (e.g. Makishima et 
al. 2001), often exhibits a significant fraction of the lumi¬ 
nosity and has to be separated from the rest of the X-ray 
emission to determine the cluster average temperature (see 
also Markevitch 199S; Arnaud & Evrard 1999| ). 

Therefore, in the ASCA data analysis, we employed a 
two-temperature (2T) model, in which isothermal plasma 
is filling the entire cluster region, while, in the central 
region, another cooler isothermal gas component is al¬ 
lowed to coexist with the hotter plasma forming a multi¬ 
phase intra-cluster medium (ICM). This 2T picture was 
estab lished for th e Ce ntaurus cluster (Fukazawa et al. 
1994 ; Ikebe et al. 1999 ) and also gave a good acco unt o f 


the ASCA data of Vi rgo/M 87 (Matsumoto et al. 
Hydra-A (Ikebe et al. 1997), Abell 1795 (Xu et al. 


1996), 


1998), 


and other nearby clusters (Fukazawa |1997; Fukazawa et 


al. pOOOl ). The hot component extends over the ma¬ 
jor part of cluster volume as far as is measurable (e.g. 
Fukazawa 1997; White 2000| ). The temperature of this hot 
component derived from the 2T model therefore represents 
the best estimation of the virial temperature of the cluster. 
Although a cooling flow spectral model also gives gener¬ 
ally a good description of the ASCA spectra (e.g. Fabian 
et al. 1994; White 2000| ; Allen et al. 2001), the estimated 
temperatures from where gas starts to cool vary consider¬ 
ably from author to author and are sometimes spuriously 
high. More importantly, recent XMM-Newton data clearly 
show that the conventional co oling f low spectral mo del is 
not adequate (Peterson et al. 2001; Tamura et al. 2001 
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Kaastra et al. 2001). The 2T model, on the other hand, 
is valid to the XMM-Newton data and the hot component 
temperature well represents the temperature of the outer 
main component (see Ikebe ^001 ). 

It is known that, apart from the central cool compo¬ 
nent, many clusters show deviations from isothermality 
such as an asymmetric temperature distribution due to 
merging (e.g. Briel & Henry 1994) or a global temper¬ 
ature decrement towards the outside (Markevitch et al. 
1998|). For such clusters, the 2T model fit gives insignifi- 


brightest region may contribute to the outer region spec¬ 
trum (Serlemitsos et al. 1995; Takahashi et al. 1995). We 


therefore, simultaneously fitted the four spectra with the 
2T model, taking into account the flux-contamination ef¬ 
fect. As a plasma code, we used the MEKAL model (Mewe 
et al. 1985, 198£; Kaastra 1992, Liedahl et al. 1995|), t hat 
is implemented in XSPEC version 10.0 (Arnaud 1996). 


cant flux to the cool component and works practically as 
an isothermal model. Thus the hot component tempera¬ 
ture gives the average temperature, which also can be a 
good measure of the virial temperature. 

Among the 106 sample clusters, 88 clusters have been 
observed with ASGA, and all the data sets have become 
publicly available by now. We retrieved the ASCA data 
from the ASCA archival data base provided by NASA 
Goddard Space Flight Center and by Leicester University. 
Here we will give a brief explanation of the analysis 
method of the ASCA data. 

ASCA has four focal plane instruments, two SIS (Solid- 
state Imaging Spectromete r) an d two GIS (Gas Imaging 
Spectrometer; Ohashi et al. 1996 ), which were usually used 
simultaneously to observe an astronomical object. In all 
the observations, the GIS were operated in the normal 
PH mode and the SIS were operated in FAINT mode or 
BRIGHT mode. The number of active CCD chips of each 
SIS used was 4, 2, or 1, corresponding to the field of view 
of 22' X 22', 22' x 11', or 11' x 11', respectively. 

We discarded the GIS and SIS data that were taken 
when the elevation angle of the X-Ray Telescope (XRT) 
from the local horizon was less than 5°. An additional 
screening requirement, that the elevation angle from the 
sunlit earth be greater than 25° and 20°, was applied to 
the SIS-0 and SIS-1 data, respectively. In order to ensure a 
low and stable particle background, we also discarded GIS 
and SIS data acquired under a geomagnetic cutoff rigidity 
smaller than 6 GV. 

The background is composed of cosmic X-ray back¬ 
ground and non-X-ray background, which were estimated 
for each cluster observation from the data of blank-sky 
observations as well as the data taken when the X-Ray 
Telescope was pointing at the dark (night) earth. 

In order to maximize the detection efficiency of the 
central cool component, we accumulated spectra over two 
regions, a central region of 2' radius and an outer region 
from 2' to Cmax, from each cluster data set. r„iax, the max¬ 
imum radius, is defined as the radius within which at least 
95% of all detected source X-ray counts, within the field of 
view of the GIS, are included. The data from the two GIS 
sensors and two SIS sensors were combined. Thus for each 
cluster, four spectra: two central and two outer spectra for 
GIS and SIS, were derived. 

An effective area as a function of X-ray energy was 
calculated for each spectrum with an ASCA simulator 
(SimARF). The point-spread-function of the ASCA XRT 
is so largely extended that the X-ray flux from the central 


A best-fit model was derived by a chi-square minimiza¬ 
tion method. The cool and hot component temperatures 
are separate free parameters. Many clusters do not require 
an additional cool component to fit the spectra, and the 
cool component temperatures cannot be constrained. For 
each cluster we investigated the significance of adding a 
cool component by comparison with a fit of a hot isother¬ 
mal model and performing an F-test. For a cluster show¬ 
ing a low significance for the additional cool component 
in the 2T model, we fixed the temperature of the cool 
component at Tiso/2, where Tjso is the temperature de¬ 
rived with the isothermal model fitting. The choice of the 
fixed cool-component temperature is based on a clear cor¬ 
relation found in the 2T model fitting that, for clusters 
showing a high significance for the cool component, the 
cool-component temperature is always very close to half 
of the hot component temperature (see Ikebe 2001). The 
hot component temperatures thus derived are summarized 
in Table 1. A more detailed analysis procedure description 
and results on the cool components will be presented in a 
following publication. 

We compared the temperatures derived here with 
those of previous worker s. The temperatures agree with re¬ 
sults of Fukazawa ( 19971 ) within 1.5% on average (Fig. |^), 
who performed a similar analysis on a smaller number of 
clusters with ASCA data, accumulating spectra from a 
central region of 2-3 arc-minutes radius and an outer re¬ 
gion, and fitting them with the 2T model. On the o ther 
hand, t emper atures given by Markevitch et al. ( 1998 ) and 
White ( 2000 ), who used different (cooling flow) models to 
fit the ASCA spectra, are higher than our values respec¬ 
tively by a factor of 1.09 and 1.25 on average for clusters 
hotter than 6 keV (Fig. |p3,c). This difference causes a sig¬ 
nificant difference in the resulting XTFs (see Sect. 3.2). 


2.3. L — T relation 

For the 88 clusters observed with ASCA, X-ray fluxes and 
luminosities were calculated from the spectral-analysis 
results described above and the ROSAT PSPC count 
rates in 50-201 PI channel given by Reiprich & Bohringer 
( ^002 ). Assuming isothermality in each cluster at the hot- 
component temperature and the Galactic column density 
given in Table 1, we determined the 0.1-2.4 keV fluxes 
and luminosities that reproduce the PSPG count rates. 
As summarized in Table 1, the luminosities were calcu¬ 
lated for two different cosmological models with an open 
(Hin,o = 0.2, Da^o = 0) and a flat (Di„_o = 0.2, Da,o = 0.8) 
universe. 
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Fig. 1. Comparison between the cluster tem peratures 
derived here and thos e by Fukazawa ( 1997 ) (a), by 
Markevitch et al. ( 1998 ) (b), and by White ( |2000 ) (c). 


We then established a correlation between the lumi¬ 
nosity in the 0.1-2.4 keV band, Lo.i- 2 . 4 keVj and the tem¬ 
perature of the hot component, T),, as illustrated in Fig. ||. 
The To.i-2.4keV"2h relations for the open and flat universe 



Fig. 2. The luminosities in the energy range of 0. 1-2.4 keV 
in the cluster rest frame in a flat universe (flm.o = 0-2, 
^A,o = 0-8) as a function of the hot component temper¬ 
ature obtained from the two-temperature model fitted to 
the ASCA data. The best-fit power-law function (eq. ||) is 
overlayed. 


were fitted individually with power-law functions by a lin¬ 
ear regression method. For the fit we did not use clusters 
with a hot component temperature cooler than 1.4 keV, 
because it is known that a single power-law function is 
not a good representation of the L — T relation over a 
wide temperature range. In particular, including clusters, 
groups, and elliptical galaxies, the L — T relat ion becomes 
steeper below ^ 1 keV (e.g. Xue & Wu 2000 ). Moreover, 
the theoretically predicted XTFs as given in Sect. 4 may 
not be directly compared with an observed XTF (see Sect. 
4.3 for more explanation). Therefore, in the power-law fit¬ 
ting, we exclude 6 sample clusters cooler than 1.4 keV 
and use only 82 clusters. For the linear regression fit, we 
employed the B CES(X 2 |Xi) estimator given by Akritas & 
Bershady ( 19961 ) (In this case, the luminosity and temper¬ 
ature are assigned as X 2 and Xi, respectively). The choice 
of the estimator is based on an argument by Isobe et al. 
( 1990| ). The best-fit functions are 

logLo.i-2.4kev(^“^ergs s“^) = 

42.31 (±0.09) ± 2.47 (±0.14) logTh(keV) (1) 


for the open universe, and 

logLo.i-2.4kev(/i“^ergs s“^) = 

42.30 (±0.09) ± 2.50 (±0.14) logTh(keV) (2) 


for the flat universe. 

We have made a correction to the derived ±o.i- 2 . 4 keV~ 
Th relation. Since the less luminous clusters are less likely 
to be sampled, the apparent mean luminosity at a given 
temperature, which is estimated by the simple regression 
fit, is biased towards higher values than the true mean 
value. As will be shown in Sect. 4, we studied this ef¬ 
fect and corrected the normalization to derive the “true” 
±o.i-2.4keV — Th relation, as 

logLo.i- 2 . 4 kev(^“^ergs s“^) = 42.15±2.47 logrh(keV)(3) 
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for the open universe, and 

logLo.i-2.4kev(h“^ergs s“^) = 42.14-1-2.50 logrh(keV)(4) 

for the flat universe. The correction reduces the lumi¬ 
nosity at a given temperature by 30%. The corrected 
Lo.i- 2 . 4 keV~ 2 h relation will also be used to calculate 14 iax 
in Sect. 3. The 30% reduction of the mean temperature 
results in a reduction of Vhiax (see Sect. 3) by 40%. 

Using the corrected To.i-2.4keV — Th relations, we esti¬ 
mated the temperatures of 14 clusters for which no ASCA 
data nor other spectroscopic measurement are available. 
An To.i-2.4keV and Th combination was chosen from the 
To.i-2.4keV ~ Th relation such that the predicted PSPC 
count rate equals the observed value. Scatter around the 
best-fit power-law in the L — T relation was assigned as er¬ 
ror. The estimated temperatures and luminosities for the 
open and flat universe cases as well as the corresponding 
0.1-2.4 keV fluxes which are the same for the open and 
flat universe cases are summarized in Table 1. 


2.4. Construction of the flux limited sample 


We set a flux-limit of 1.99 x 10“^^ ergs s“^ cm“^ in 0.1- 
2.4 keV band to construct a flux-limited complete sample 
from the master catalog of 106 clusters. As in Reiprich 
& Bohringer ( ^002 ) for HIFLUGCS, additional selection 
criteria applied are that the absolute Galactic latitude is 
greater than 20°, and the cluster is located neither in the 
Magellanic Clouds nor in the Virgo cluster region. The 
total number of clusters thus selected is 63, whose tem¬ 
peratures and redshifts are shown in Fig. Since the 
temperatures of the individual clusters were newly deter¬ 
mined in this paper, the fluxes estimated from the PSPC 
count rate ar e sligh tly revised from the value in Reiprich 
& Bohringer ( 2002| ) where the temperatures are compiled 
from literature. All of the 63 clusters selected here, how¬ 
ever, still correspond with the 63 HIFLU GCS clusters. 
Tests performed in Reiprich & Bohringer ( ^002 ) showed 
no indication of incompleteness of the sample. 

In the following analysis, we exclude 2 clusters cooler 
than 1.4 keV (see Sect. 2.3 and 4.3) from the flux-limited 
complete sample, and use 61 clusters to study the XTF 
and make constraints on cosmological parameters. This is 
the largest complete sample of clusters up to now with 
temperature measurements or reliable temperature esti¬ 
mates. We also stress the accuracy of the temperature 
measurements for the sample clusters: Among the 61 clus¬ 
ters used for constructing the XTF, the temperatures of 
56 clusters (90% of the sample) were measured with ASCA 
data with 90% errors of 6-10%, other spectroscopic mea¬ 
surements were used for 3 clusters, and there are only two 
clusters for which the temperature is estimated with the 
L — T relation. 


3. The Temperature Function 



Fig. 3. X-ray temperature vs redshift for the clusters of 
the flux-limited sample. 


3.1. Derivation of XTF 


Using the X-ray flux-limited complete sample of 61 clus¬ 
ters (Th > 1.4), we address in this section the XTF de¬ 
scribing the number density of clusters of galaxies in co¬ 
moving coordinates at present as a function of the X-ray 
temperature. A full description of the analysis method 
can be found in the Appendix. However, the key points 
are presented here. The differential temperature function, 
4>{T), is defined as the number density of clusters in the 
temperature range from T — ^ to T + ^. From a finite 
number of sample clusters, it can be evaluated as 


ct>{T)dT = 

T-.^<Ti<T+^ 



(5) 


where i denotes a sequential number of clusters in a tem¬ 
perature range from T — ^ to T -|- and Umax(Ti) is 
a maximum comoving volume in which a cluster having a 
temperature of T could have been detected under the se - 
lection criteria given above. Following Markevitch ( 1998 ), 
we assume, for a given temperature, the logarithmic lumi¬ 
nosity, log T, follows a Gaussian distribution around the 
mean which is given with the L — T relation, i.e. T = AT°^, 
and with a constant standard deviation of criogi- Umax(T) 
then can be evaluated as 




(T) = 

W(T,T) 


27rCTiogL 


exp • 


(logAT“-logT)2] ^ 

- ^ - '> dlogT, (6) 




where Uniax(T,T) is the maximum search volume for a 
cluster having the luminosity of T and the temperature of 
T and is given by Eq. A.9. Recall that 14iax(T) depends 
on the cosmological parameters and was calculated for an 
open (Um^o = 0.2, Ua,o = 0) and a flat (Um^o = 0.2, LIa,o = 
0.8) universe. The results are listed for individual clusters 
in Table 2. The XTF thus derived with 1 keV bin widths 
is illustrated in Fig. 

The improvement of the new result over previous work 
is mainly characterized by a significant broadening of 























6 


Ikebe et al.: A New Measurement of the X-ray Temperature Function of Clusters of Galaxies 



Fig. 4. The X-ray temperature function derived with the 
flux-limited complete sample of 61 clusters for the case 
of a flat universe = 0.2,nA,o = 0.8). Filled circles 

show the XTF evaluated with 1/Fmax(7i), while open stars 
show the one obtained with l/uniax(Fi) 21). The tempera¬ 
ture bin widths are 1 keV, except for the lowest tempera¬ 
ture bin that includes fcT=1.4-2.5 keV. The vertical error 
bars for filled circles indicate Poisson errors. Open stars 
should have the same Poisson errors, which are omitted 
for a clear display. The solid curve is the best-fit Press- 
Schechter function derived in Sect. 4.1. 

the temperature range towards lower temperature sys¬ 
tems down to 1.4 keV. A single power law is no longer 
a good description to the observationally measured XTF 
and an exponential cut off towards the higher tempera¬ 
ture is clearly seen, as a theoretical model predicts. The 
detection of such a curvature brings tight constraints on 
cosmological parameters as will be shown in the next sec¬ 
tion. 

Conventionally, 14iax(2i) in Eq. ^ is often simply re¬ 
placed by Uniax(Fi,Ti) to evaluate the XTF (e.g. Henry 
|2000| ). For comparison, the XTF evaluated with (f>{T)dT = 

l/'*^max(Fi, 2i) is also overlayed in Fig. H, which overall 
shows good agreement with the one with X]l/^ax(2i)) 
but exhibits larger variance due to the intrinsic scatter in 
the L — T relation. 


3.2. Comparison to previous results 

In Fig. we compare our result to previously obtained 
XTFs overlayed in the form of best-fit power law func¬ 
tions. In the 3 — 8 keV range, our XTF shows good agree¬ 
ment with the previous results. 

We have derived a cumulative XTF, which is the num¬ 
ber density of clusters that are hotter than a certain tem¬ 
perature, by 


n{T) = 

T<Ti 


1 


l^max(2i) ■ 


(7) 


Table 2. Summary of the flux-limited complete sample 


NAME 

fcrh(keV) 

Open 

Mpc-3) 

Flat 

A2163 

10.55 

2.39e+08 

2.65e-b08 

A0754 

9.00 

1.40e+08 

1.52e-f08 

A2142 

8.46 

1.14e+08 

1.23e-f08 

COMA 

8.07 

9.69e-b07 

1.04e-f08 

A2029 

7.93 

9.12e+07 

9.78e-f07 

A3266 

7.72 

8.30e+07 

8.89e-f07 

A0401 

7.19 

6.50e+07 

6.91e-f07 

A0478 

6.91 

5.66e+07 

6.00e-f07 

A2256 

6.83 

5.45e+07 

5.77e-f07 

A3571 

6.80 

5.36e+07 

5.67e-f07 

A0085 

6.51 

4.60e+07 

4.87e-f07 

A0399 

6.46 

4.46e+07 

4.72e-f07 

A2204 

6.38 

4.28e+07 

4.52e-f07 

ZwC11215 

^ 6.30 “ 6.36 

4.09e-b07 

4.49e-f07 

A3667 

6.28 

4.05e+07 

4.28e-f07 

A1651 

6.22 

3.92e+07 

4.14e-f07 

A1795 

6.17 

3.81e-b07 

4.02e-f07 

A2255 

5.92 

3.29e+07 

3.46e-f07 

A3391 

5.89 

3.23e+07 

3.40e-f07 

A2244 

5.77 

3.01e-b07 

3.17e-f07 

A0119 

5.69 

2.86e+07 

3.00e-f07 

A1650 

5.68 

2.85e+07 

2.99e-f07 

A3395S 

5.55 

2.64e-b07 

2.76e-f07 

A3158 

5.41 

2.40e+07 

2.50e-f07 

A2065 

5.37 

2.34e+07 

2.44e-f07 

A3558 

5.37 

2.35e+07 

2.45e-f07 

A3112 

4.72 

1.48e+07 

1.53e-f07 

A1644 

^ 4.70 

1.46e-b07 

1.51e-f07 

A0496 

4.59 

1.34e+07 

1.39e-f07 

A3562 

4.47 

1.22e-b07 

1.27e-f07 

A3376 

4.43 

1.18e+07 

1.22e-f07 

A2147 

4.34 

l.lOe+07 

1.13e-f07 

A2199 

4.28 

1.04e+07 

1.08e-f07 

A2597 

4.20 

9.75e+06 

1.01e-f07 

A0133 

3.97 

8.02e+06 

8.25e-f06 

A4059 

3.94 

7.76e+06 

7.98e-f06 

A0576 

3.83 

7.07e-b06 

7.26e-f06 

HYDRA-A 

3.82 

6.98e+06 

7.16e-b06 

A3526 

3.69 

6.15e+06 

6.30e-|-06 

A1736 

3.68 

6.11e-b06 

6.26e-|-06 

2A0335 

3.64 

5.87e-b06 

6.00e-b06 

A2063 

3.56 

5.41e-b06 

5.53e-f06 

A1367 

3.55 

5.36e+06 

5.48e-f06 

A2657 

3.53 

5.25e+06 

5.36e-f06 

A2634 

3.45 

4.84e-l-06 

4.94e-|-06 

MKW3S 

3.45 

4.87e-b06 

4.97e-f06 

A2589 

3.38 

4.48e+06 

4.57e-f06 

MKW8 

3.29 

4.08e-b06 

4.16e-f06 

A4038 

3.22 

3.78e+06 

3.85e-f06 

A1060 

3.15 

3.50e+06 

3.56e-f06 

A2052 

3.12 

3.37e-b06 

3.43e-f06 

IIIZw54 

^ 2.98 * 3.00 

2.87e-b06 

2.99e-f06 

EXO0422 

^ 2.90 

2.60e-|-06 

2.63e-|-06 

SllOl 

^ 2.60 

1.76e-b06 

1.77e-f06 

A0400 

2.43 

1.37e-b06 

1.38e-f06 


The resulting cumulative XTF is illustrated in F ig. || 
together with other recent results by Henry (200C|) and 
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Table 2. continued 


NAME 

fcrh(keV) 

Open 

Mpc-3) 

Flat 

A0262 

2.25 

1.04e-t06 

1.04e-t06 

MKW4 

1.84 

5.00e-t05 

4.96e-t05 

A3581 

1.83 

4.88e-t05 

4.84e-t05 

FORNAX 

1.56 

2.72e-t05 

2.68e-t05 

N1550 

1.44 

2.05e-t05 

2.02e-t05 

N507 

1.40 

1.83e-t05 

1.79e-t05 

N5044 

1.22 

1.10e-t05 

1.08e-t05 

N4636 

0.66 

1.16e-t04 

1.10e-t04 


: as Table 1. 


Markevitch (1998). Our result shows agreement with 
Henry’s XTF constructed from his sample of 25 clusters. 
The individual temperatures of the 25 clusters were eval¬ 
uated by averaging the previous temperature measure¬ 
ments with ASCA and other missions based on single¬ 
temperature fits (see Henry 200C), which agree with the 
hot component temperatures obtained in this work on av¬ 
erage within 1%. Above 6 keV, there is apparent disagree¬ 
ment between our result and the result by Markevitch 
(1998) showing a higher amplitude. This is due to a larger 
number of clusters in the highest temperature bins in the 
Markevitch sample, despite the fact that smaller t4iax val¬ 
ues apply to his sample. (Markevitch selected his sample 
with the same flux limit as ours but also limited the red- 
shift range to 0.04 < z < 0. 09). T he larger number of very 
hot clusters in Markevitch ( 1998 ) is due to higher temper¬ 
ature estimates compared to our values for temperatures 
above ~ 6 keV (Fig. ^). 


Blanchard et al. (EOOO), using a sample and tempera¬ 
tures mainly from the XBACS (Ebeling et al. 1996D and 
the ASCA measurement by Markevitch et al. (1998), re¬ 


spectively, obtained the XTF, which virtually perfectly 
agrees with that by Markevitch ( 199^ . Pierpaoli et al. 
(2001), who have made another attempt to derive the 
XTF, used the same sample a s Mar kevitch (199S) but the 
temperatures given by White (|2000| ). The Pierpaoli cumu¬ 
lative XTF has the highest amplitude among other works. 
This is bec ause of the systematically higher temperatures 
by White (200C) (see Fig. ||). 


4. Constraint on cosmological parameters 

In the following, using the flux-limited complete sample 
of 61 clusters (Th > 1-4) and theoretical XTF models, we 
obtain constraints on cosmological parameters. The XTF 
at the present epoch is most sensitive to the amplitude of 
density fluctuations on a scale of 8/i“^ Mpc, as, and least 
sensitive to the cosmological constant, A. We therefore 
concentrate on two different Ha,o model families includ¬ 
ing open (flm_o < 1, = 0) and flat (Hm^o + f^A^o = 1) 

universes. Recent measurements of high-redshift SN la in¬ 
dicate a positive cosmological constant (Perlmutter et al. 



Fig. 5. The X-ray temperature function derived with the 
flux-limited complete sample of 61 clusters for the case of 
a flat universe (Hm.o = 0.2, Ha, o = 0.8). The three lines 
indicate previous measurements in the form of the best- 
fit power-law functions. Solid, dotted, a nd da shed lines 
corr espon d to the resu lts by Markevitch ( 1998 ), Edge et 
al. ( |l990 ), and Henry ( ^000[) , respectively. 



Fig. 6. The cumulative temperature function for the case 
of a flat universe (Hm.o = 0.2, Ha,o = 0.8) is shown in 
grey solid line with its 68% error band. The errors are es¬ 


timated by ■\J'^Ti<T (l/i^ax(’71))^. Previous results from 
Markevitch ( 1998| ) and Henry ( pOOO ) are also shown by 
the solid and dotted black lines, respectively. 


1999; Riess et al. |1998 ). Moreover, the cosmic microwave 
background (CMB) temperature-fluctuation map recently 
derived with the Boomerang experiment stron gly supports 
a flat universe model (de Bernardis et al. 2000 ). Therefore, 
a flat universe model with a positive cosmological constant 
is currently the most plausible description of our universe. 
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For the purpose of comparison with previous works, we 
also consider the case of an open universe. 

For the theoretical prediction of the XTF, we employ 
two different models, a conventional Press-Schechter mass 
function (Sect. 4.1) and the Formation-Epoch model by 
Kitayama & Suto ( 1996| ) that is a modification of the 
former theory (Sect. 4.2). The model XTF depends also 
on the Hubble constant, which is here assumed to be 71 
km/s/Mpc (i.e. h = 0.71) based on the latest results of 


given in Kitayama & Suto ( |1996 ). Since Rv/Rc >> 1, the 
virial mass is approximately given as, 


= ^(ttpoA) 


P 


firUpG^ 


3/2 


kT 


3/2 


( 12 ) 


the HST Hubble constant project (see Mould et al. 2000 
and references therein). 

4.1. Press-Schechter model 

According to the Press-Schechter formalism, the number 
density of clusters with a mass in the range of M to M -|- 
dM that have collapsed before a redshift z is given as 


2 ps{M, z)dM = 
5c{z) 

TT M 


da{M) 


dM 


exp - 


Sliz) 

2 cr 2 (M) 


dM 


( 8 ) 


where po is the mean density of the universe at present, 6c 
is the linear overdensity of a system collapsed at a redshift 
z and evaluated at present (for which we use the formula 
given in Kitayama & Suto 19961 ), and cr^ is the variance of 
mass extracted from the density fluctuation field within 
a spherical region that includes a mass M on average. 
For the mass variance, we empl oyed an analytical formula 
derived by Kitayama & Suto ( 1996 ) based on the CDM 
power spectrum, which is given as 


Note that when P/Pspec = 2/3, {Pspec = ), this 

equation becomes identical to that derived from a mod¬ 
ified singular isothermal sphere by Bryan & Norman 
(1998). Although an additional modification to the M — T 
relation is often applied by multiplying a factor accounting 
for departures from virial equilibrium (e.g. Henry 2000), 
we did not introduce this factor here. The observed values 
of P are distributed in a relatively narrow range around 
~ 2/3 and may show a weakly positive correlation with 
temperatures (Arnaud & Evrard 1999; Vikhlinin et al. 
|1999| ; Horner et al. 1999). Vikhlinin et al. measured P 
values systematically with iJOS'AT'PSPC data excluding 
a central region inside 30% of the virial radius and found 
a correlation of P with X-ray temperature, which can be 
fitted with the linear function /3 = 0.017 T(keV) -I- 0.621 
(Henry 2000). We employ this relation and substitute it 
into Eq. |12|. The effect of using different M — T relations 
will be studied in Sect. 5. 

We further assume that a cluster collapsed just be¬ 
fore it is observed. This recent formation approximation, 
which has been extensively applied in previous conven¬ 
tional modelings, is relaxed in the Formation-Epoch model 
described in the next section. Combining all the formulae 
given above, the model XTF at a given redshift is given 

by 

dMv 


(T oc (1 + 2.208 toP - 0.7669to2p + 0.7949m3P)2/(9p), (9) z) dT = npPM^{T, z), z) ^(T, z) dT . (13) 


where p = 0.0873, m = M(Th)^//10^^. T is the 
shape parameter and is given as 


Hereafter we refer to this model XTF as the PS model. 


r = tl^^oh{To/2.7K) ^exp[-OB(l + (10) 4.2. Formation-Epoch model 


where Tq and Hb are the temperature of the cosmic mi¬ 
crowave background and the baryon density, which are 
assumed to be 2.726 K (M ather et al. 1994 ) and Hb = 
0.0193/i“^ (Buries & Tytler 1998 ), respectively. 

The mass function given above is converted to the XTF 
by means of a mass-temperature relation. We assume that 
the ICM in a cluster i s desc ribed by an isothermal [3 model 
(e.g. Forman & Jones 1982 ), in which the density profile is 
given as Ugas = ngas,o[l + (i?/i?c)^]~^^^^, and that hydro¬ 
static equilibrium is achieved. The virial mass is defined 
as 


Mv = --nRlpoil -f z)^A = 


3pkTR^ {R^/Rcf 

prUpG l + (i?v/i?c) 


2 ’ 


( 11 ) 


In the model XTF given in the previous section, the for¬ 
mation redshift, zy, of an observed cluster is assumed to 
be the same as the observed redshift (i.e. recent forma¬ 
tion approximation). However, this assumption is not rea¬ 
sonable in particular in a low density universe. Kitayama 
& Suto ( 19961 ) constructed a model XTF taking into ac¬ 
count the distribution of the dark-matter-halo formation 
epoch as well as the evolution of temperature (hereafter 
FE model). They showed that the XTF predicted by the 
FE model can be significantly different from the PS model 
described in Sect. 4.1, in particular when the temperature 
evolves after the collapse. 

The FE model is given as. 


where T is the temperature of ICM, i?v is the virial radius, 
Rc is the core radius, p is the mean molecular weight of 
the ICM (set to 0.6), nip is the proton mass, G is the 
gravitational constant, and A is the ratio between a cluster 
mean density and the mean density of the universe at the 
cluster-formation redshift z. For A we use the formula 


^OO 

0FE(^obsi Z^dTcifg = / dZf 

J z 

nps(2Mv(T, zf), z)' 


dv 

2^(zf,2Mv(r,Zf),z) 






dTobs , (14) 


K(zf,z) dT 

where rips and Mv are given in Eqs. ^ and |l^ , respectively, 
and dp/dzi given by Lacey & Cole ( 1993 ) is a differential 
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distribution function of halo formation epochs of a cluster 
that has a mass of 2Mv at the observed redshift z. At the 
collapse redshift Zf, the cluster is assumed to be formed 
with a mass of M^, increasing in mass after its formation. 
K(zf, z) = gives the evolution in temperature via 

Tohsiz{,z) = K{zf, z)r(zf), where T(zf) is the virial tem¬ 
perature at the collapse redshift Zf, while Tobs(-2f, z) is the 
observed temperature at the redshift z. s = 0 corresponds 
to no temperature evolution, while s > 0 and s < 0 in¬ 
dicate positive and negative temperature evolution after 
the collapse, respectively. 


4.3. Fitting method 

Instead of fitting the model XTF to the XTF derived in 
Sect. 3, we use the predicted XTF and the L — T relation 
to calculate the expected luminosity and temperature dis¬ 
tribution for the given survey selection function which is 
given in a form of Vmax, and perform a fit to the observed 
luminosity-temperature distribution of the 61 clusters se¬ 
lected in Sect. 2.4. In this analysis, the intrinsic luminos¬ 
ity distribution for a given temperature, which in other 
words is the conditional luminosity function at a given 
temperature, is assumed to follow a Gaussian function in 
logarithmic scale with a constant standard deviation and 
a mean luminosity given by a power-law function of the 
temperature. This assumption can be justified by obser¬ 
vational results as shown in Sect. 3.1. In order to simplify 
a model calculation, we assume that the XTF does not 
change within the search volume of our flux limit and use 
a model XTF at the median redshift of our cluster sample 
at z=0.046 as a representative value. Systematic errors on 
the final results caused from this assumption are found 
to be negligibly small compared with statistical errors. 
Incorporated by another Gaussian that represents tem¬ 
perature measurement errors, the expected cluster num¬ 
ber density in a unit logarithmic temperature and a unit 
logarithmic luminosity is given as 


N{L,T) d\ogL dlogT = 


/ OO 

dr 

-CX5 


— OO 

1 


(j){T',z = 0.046) 


{L,r) 




: exp 


{h = 0.71)3 
^ (log AT°‘ — log L)^ 






: exp 


r (r-T)^] 

dT 

1 24 J 

dlogT 


dlogL dlogT, (15) 


A and a specifying the power-law function of the L — T 
relation that gives the mean of the luminosity; and diog l 
that is the constant standard deviation of the luminosity 
from the L — T relation in logarithmic scale. Although 
'I'max varies with the assumed model universe, i.e. with 
the values of flm.o as well as erg which are the fitting pa¬ 
rameters, we used a specific Umax value calculated with 
(flm.O: ^^A,o) = (0.2, 0), or (0.2,0.8), in the case of an open 
or flat universe, respectively. Since all the clusters in our 
sample are nearby, Umax is virtually constant for different 
f^m.o values, using the fixed Umax does not introduce any 
considerable systematic error on the final constraints on 
flm.o and flA^Q. The model cluster distribution was then 
fitted to the observed temperature and luminosity distri¬ 
bution. Assuming a Poisson distribution of cluster number 
counts at a given luminosity and temperature, we defined 
a logarithmic likelihood function as 

In/: = ^ln7V(ri,Li) - j N{T, L) dlogL dlogT , (16) 


where i denotes each sample cluster and the integration is 
performed over the logL(ergs s“^) range from 41.5 to 45.5 
and the logT(keV) range from 0.14 to 1.05 (corresponding 
to T = 1.4 - 11.2 keV). By maximizing the likelihood 
function, we derived the best-fit model XTF as well as 
the L — T relation. 

There are two reasons for setting the low temperature 
cut-off at 1.4 keV. First, as discussed in Sect. 2.3, the L—T 
relation breaks down for elliptical galaxies which have a 
lower gas mass fraction, and actually g ets ste eper below 
^ 1 keV as shown by e.g. Xue & Wu ( 200C| ). Thus the 
simple power-law modeling for the L — T relation is not 
appropriate for the whole temperature range. The second 
reason is based on the more essential problem concerned 
with the comparison between the theoretical XTFs and 
the observation. The systems at the low temperature end 
are mostly single elliptical galaxies and some of them may 
not be isolated but surrounded by larger scale dark matter 
halos. As has been shown by e.g. Matsushita (1998), X-ray 
brightness profiles of elliptical galaxies have large variety 
suggesting a significant variation in total mass for a given 
temperature, i.e. a considerable scatter in the M — T rela¬ 
tion would be expected. Therefore, a cluster abundance at 
the low temperature end probably does not directly corre¬ 
spond to an abundance of dark matter halos at a certain 
mass scale. 


4.4. Results 


where (j) is calculated for h = 0.71 and normalized with 
the /i“3 factor to a case of the h = 1 universe. 'Cmax(7v, T) 
is a maximum comoving volume calculated for = 1, 
in which a cluster having a temperature of T and a lu¬ 
minosity of L could have been detected under our sam¬ 
ple selection crite ria, for which the calculation method is 
given by Eq. . The temperature measurement errors 
are ctt = 0.051 -|- 0.0014 r^.sss q^j- ascA results. 
Free parameters to a fit are; Ug and Llm.o specifying (j); 


We first present the results from the PS model fitting. The 
best-fit values and errors of the fitting parameters are sum¬ 
marized in Table 3. The best-fit distribution function in 
the logL—logT space in the flat universe case is illustrated 
in Fig. |a The corresponding model XTF is also illustrated 
in Fig. y. The constraint in flm.o ~ <^8 space is shown in 
Fig. g for each open and flat universe case. The contours 
represent the 90% confidence range for two parameters of 
interest, and the constraint is notably tight. 
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In order to demonstrate an improvement on the con¬ 
straint by the widest temperature range and the largest 
sample ever achieved, we also performed the fitting with 
the PS model only using 51 sample clusters above 3 keV. 
As shown in Fig. the improvement is significant, in par¬ 
ticular, for flm.o- Since flm.o is tightly related to the over¬ 
all shape of the mass variance, i.e. the power spectrum, a 
larger mass range to measure the mass fluctuation ampli¬ 
tude is essential. 


The “true” L — T relation employed in Sect. 2.3 and 
Sect. 3.1 was estimated by this fitting. Fixing the a value 
as obtained in Sect. 2.3, to 2.47 and 2.50 for the open 
and flat universe case, respectively, we performed the PS 
model fitting and derived logA=42.15, a = 0.24 for the 
open universe case, and log A=42.14, cr = 0.24 for the flat 
universe case, which are given in Eqs. ^ and ^ and are 
applied in Eq. ^ 

The FE model gives somewhat different results. In 
addition to f2m,o and cts, the FE model introduces an¬ 
other parameter, s, that models a temperature evolution, 
which can be a free parameter in a fit. However, follow¬ 
ing Kitayama & Suto ( 1996 ), in the following analysis, we 
rather fixed s at s = 0 and s = 1 such that they represent 
no evolution and a maximum possible positive evolution 
of the temperature, respectively. Since there is no obser¬ 
vational evidence for the evolution of the temperature, 
e.g. no redshift dependency in the L — T relation is found 
(Mushotzky & Scharf 1997; Donahue et al. 1999; Della 
Ceca et al. ^000 ), the no evolution model may be more 
plausible. However, the effect of the temperature evolu¬ 
tion is still interesting to study here. The fitting results 
are summarized in Table 3, and the constraints on Dni,o 
and tTg are illustrated in Fig. The FE model in the no 
temperature evolution case (s = 0) shows similar results 
with those of the PS model, since the XTF predicted by 
the s = 0 FE model is very clo se to that of the PS model 
(see Pig.8 in Kitayama & Suto 1996 ). On the other hand, 
in the maximum temperature evolution cases (s = 1), the 
FE model gives significantly smaller erg values than those 
with the PS model. 


5. Discussion 


5.1. Uncertainty in the M — T relation 

The M — T relation used in our analysis may be sub¬ 
ject to uncertainties introducing a significant systematic 
error in the final results on the cosmological parame¬ 
ters. Theoretical argu ment s yield a scaling law, that is 
M oc r3/2 Kaiser 1986| ). A number of numerical sim¬ 
ulations support this scaling law, but different normaliza¬ 
tions are derived from different simulations. Henry ( ^000 ) 
summarized such differences in the normalization of the 
M—T relation derived from different hydrodynamical sim¬ 
ulations. A cluster mass may vary from 62% to 119% of 
the mean value. 

The M — T relations derived observationally from X- 
ray mass measurements give somewhat different results. 



log Th(keV) 


Fig. 7. The best-fit model function in the luminosity- 
temperature space for the case of the flat universe is shown 
with contours. The 63 clusters from the flux-limited com¬ 
plete sample are overlayed with open circles. The vertical 
dotted lines indicate the temperature range used for the 
fitting. Therefore, 61 clusters are used to make constraints. 


Horner et al. ( 199£ ) derived an M — T relation which has 
a similar exponent as the scaling law but a lower normal¬ 
ization than the numerical simulations. On the other hand, 
Nevalainen et al. ( ^000 ) and Finoguenov et al. ( 2001 ) ob¬ 
tained steeper relations of and respectively, 

and lower normalizations. 

In order to investigate the effects of the different M — T 
relations, we repeated the fitting performed in Sect. 4 with 
the PS model for the flat universe case. The six M — T re¬ 
lations we have examined are summarized in Table 4. Each 
M — T relation can be expressed with Eq. |^, replacing 
the [3 value with the formula given in Table 4. In Fig. ^ 
the best-fit parameters derived with the different M — T 
relations are compared with the result presented in Sect. 
4.4. The effect of using different M — T relations is signifi¬ 
cant and the systematic errors introduced are comparable 
to, or even larger than, the current statistical errors. As 
seen in Fig. ^ a change in the index of the M — T re¬ 
lation moves the best-fit point along the elongation direc¬ 
tion of the confidence contour, while an amplitude change 
moves the best-fit point along the perpendicular direction. 
Taking into account these systematics from the different 
M—T relations, the constraints on f2ni,o and erg are revised 
and summarized in Table 5. 


With the FE model, where the recent formation ap¬ 
proximation is not valid, it is not appropriate to use the 
empirical M — T relation obtained based on the observa¬ 
tions of nearby clusters. This is because the current M — T 
relation should be a result of accumulation of different 
M — T relations from different formation redshifts which 
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Fig. 8. Constraints on flmfi and as with the PS model 
in case of the open (upper panel) and flat (lower panel) 
universe. Bold contours are results including all sample 
clusters hotter than 1.4 keV, while dotted contours are 
derived when only clusters hotter than 3.0 keV are used. 

Table 4. M-T relations 


Source 

P 

Hydrodynamic simulation (mean) 

1 X 1.21 

Hydrodynamic simulation (maximum) 

1 X 1.36 

Hydrodynamic simulation (minimum) 

1 X 0.88 

Horner et al. 

0.496 

Nevalainen et al. 

0.437 

Finoguenov et al. 

0.415 


depend on the cluster masses. Therefore, the M — T re¬ 
lations obtained from the Hydrodynamic simulations may 
be more appropriately applied with the FE model. 


Fig. 9. Constraints on O„i,o and as in case of the open 
(upper panel) and flat (lower panel) universe. Thin closed 
contours show results with the PS model. Bold dotted and 
solid closed contours are derived with the s = 1 and s = 
0 FE models, with and without a temperature evolution 
effect, respectively. The line shows the co nstrain ts from 
the CO BE normalization by Bunn & While ( 1997 ) for h — 
0.71. 


Table 5. Results with different M-T relations. The pa¬ 
rameter ranges correspond to 90% statistical error for one 
parameter of interest plus systematic errors with different 
M — T relations applied. 


Model XTF 


,0 

as, 


PS(Open) 

0.09 - 

0.37 

0.55 - 

1.13 

FE(Open,s=0) 

0.07 - 

0.32 

0.61 - 

1.15 

FE(Open,s=l) 

0.03 - 

0.23 

0.52 - 

0.80 

PS(Flat) 

0.10 - 

0.38 

0.57 - 

1.19 

FE(Flat,s=0) 

0.10 - 

0.35 

0.65 - 

1.33 

FE(Flat,s=l) 

0.06 - 

0.18 

0.66 - 

0.97 


5.2. Comparison with other measurements 


We compared our constraints on th e cos mological param¬ 
eters with previous works. Bahcall (200C) summarized the 
cluster constraints on flm.o and as as 0.2±0.1 and 1.2±0.2 
(68% confidence level), respectively. These values are con¬ 
sistent with our results. 


The power spectrum of matter distribution in the 
present universe has been measured from the spatial dis¬ 
tribution of galaxies and galaxy clusters, and found to 
be generally well repro duced by the CDM pow er spec¬ 
trum (Peacock & West 1992; Einasto et al. 1995; Jing & 
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Fig. 10. The best-fit values of flm,o and as derived with 
the PS model in the case of the flat universe, for which 
different M — T relations are applied. A 90% conhdence 
contour derived in Sect. 4 is also overlayed. 


Valadarnini I993|; Einasto et al. 1997; Retzlaff et al. 1998; 
Tadros et al. 1998 ; Miller & Batuski 2001 ). An important 
feature of the power spectra is a turnover, where the scale 
is closely linked to the horizon scale at matter -radiation 
equality. Most recently Schuecker et al. ( 2001 ), using a 
sample of 452 X-ray selected clusters of galaxies, deter¬ 
mined the location of the maximum of the power spectrum 
in the range of 0.022 < /cmax < 0.030 h Mpc“^. According 
to the CDM power spectrum, fcmax is well approximated 

by 


fcmax = 0.114 r h MpC 


-1 


(17) 


where P, given in Eq. |^, is a function of flm.o- Therefore 
one finds 0.33 < flm.o < 0.44. The constraints on flm.o we 
have derived in Sect. 5.1 with the PS model yield 0.005 < 
fcmax < 0.027 h Mpc“^, which overlaps with the result 
by Schuecker et al. Our results with the FE model give 
4 X 10“"^ < fcniax < 0.022 h Mpc“^ in the case of the open 
universe, and 2 x 10“^ < k-^sx < 0.025 h Mpc“^ in the 
case of the flat universe, both of which are also consistent 
with the result by Schuecker et al. 

A measurement of the matter density fluctuation at 
the largest spatial scale was performed using the COBE 
DMR experiment (e.g. Bennett et al. 1996). Bunn & 
White ( 1997 ), using the CDM power spectrum model, ex¬ 
trapolated the COBE results to cluster scales and gave 
constraints in Dm.o-o’s space, which are shown in Fig. 

In the case of the open universe, our results with any of 
the model XTFs do not overlap the COBE constraints, 
even if the effects of the uncertainty in the M — T relation 
are taken into account. On the other hand, in the case of 
the flat universe, our results with the effects of the un¬ 
certainty in the M — T relation, are consistent with the 
COBE constraints regardless of the choice of the model 
XTF. 


6. Summary and Conclusions 

The most important result concerning the XTF derived 
here is that it covers a large enough range of tempera¬ 
tures to reveal the shape of the function beyond a simple 
power law representation. As expected from theoretical 
considerations, there should be a very sharp cut-off at 
the high mass end and the function should turn into a 
shallower slope at the low mass end as shown by our ob¬ 
servations. We have demonstrated that these constraints 
on the actual slope of the XTF allow us to derive inde¬ 
pendent constraints on the two parameters which are most 
important in cosmic structure formation models: and 

as- Leaving out only 16% of the sample clusters from the 
low temperature end almost doubles the uncertainty for 
these independent constraints. Therefore, while these in¬ 
dependent constraints are an important achievement of 
the present work on one hand, this demonstration also 
shows that these constraints are sensitive to any system¬ 
atic error introduced into the XTF. 

For this reason we have studied the stability of our 
results to variations introduced to the observed data and 
also variations of the theoretical modeling. We have fully 
incorporated the measurement errors in fluxes and tem¬ 
peratures of the sample clusters. With the typical errors 
of 5% for the flux and 6-10% for the temperature, the bi¬ 
ases were found to be relatively small and not to cause 
a drastic change on the hnal results in our analysis. On 
the other hand, we have shown that there are considerable 
variations and uncertainties in the theoretical modeling of 
the XTF, which introduce significant systematic errors in 
the hnal constraints on Q^.o and as- Various M — T rela¬ 
tions that are employed to model the XTF have resulted 
in considerably different values of Oi„,o and as- Another 
source of variation in the theoretical modeling that we 
have studied is a recent formation approximation which 
is employed with the Press-Schechter formalism to build 
the PS model here. Using the FE model, in which the dis¬ 
tribution of cluster formation redshifts and the evolution 
of the mass as well as temperature of a cluster after the 
collapse are analytically formulated, we have shown that a 
considerable difference could be introduced in the derived 
cosmological parameters, in particular in as- This further 
motivates us to establish a more sophisticated model to 
the XTF incorporating new X-ray observations of the high 
redshift clusters with Chandra and XMM-Newton. 

Finally, taking into account all the uncertainties, we 
put our general constraints on and as as f2ni,o=0.03- 
0.37, (78=0.52-1.15 in an open universe and Uni,o=0.06- 
0.38, (78=0.57-1.33 in a flat universe, respectively. From 
the comparison with the COBE constraints, a flat universe 
is more preferable than an open universe. 


Appendix A: Calculation of Umax 

Here we present our method, used in Sect. 3, to evalu¬ 
ate the X-ray temperature function in more detail. A flux 
limited sample has been extensively used for constructing 
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the X-ray luminosity function (XLF) by many authors 



present the first description of how to construct this func¬ 
tion in the presence of measurement uncertainties, correla¬ 
tion scatter and selection biases. For a sample with a finite 
number of objects, the differential XLF, i.e. the number 
density of clusters having the rest frame luminosity in a 
range of L — ^ to L + ^ defined in a comoving space, 
can be evaluated as 


(/)(L)AL = 


E 


L-i 


-<Li<L + i 




(A.l) 


where i denotes individual clusters in the luminosity range 
of a bin. v„iayi{L,T) is a maximum search volume mea¬ 
sured in comoving space where a cluster having the rest 
frame luminosity of L and the temperature of T could 
have been detected under the flux limited condition. 
Introducing ^max, the maximum redshift at which such a 
cluster could have been detected, and I4om, a total volume 
in the comoving coordinate from us to a given redshift, we 
can derive Vm&x{L,T) as 

^max(.b,T') = - Vj;om ('2'max) j (-^■^) 

47r 


where is the total sky coverage in steradians, and Va 
is, in the case of a flat universe, simply given by 

3 




47r f Di{z) 
1 + z 


(A.3) 


where Di is a luminosity distance, while, in the case of an 
open universe, it is given with a more complicated form 
(see e.g. Sandage 19881 ). z^ax is given by solving the fol¬ 
lowing equation 


.^l(^max) — 


Lohs{L, r, •^^max ) 
47r/obs ,lim. 


(A.4) 


where /obs.iim is the flux limit defined in a certain en¬ 
ergy range in the observed frame, and Lobs is the same 
observed-frame luminosity that is expected if the cluster 
is located at Zmax- In order to solve this equation with 
the LT-correction, the spectral energy distribution, i.e. the 
temperature, T, is required, although the dependency of 
-^max on T is generally weak. 

Using u„iax(L, T), the differential XTF can be given as 


(j){T)AT = 


‘^max(Li, T)) 


(A.5) 


This estimator is used in e.g. Henry (2000). However, since 
there is a distribution of luminosities for a given temper¬ 
ature, this method introduces a scatter in a derived XTF, 
unless there are large numbers of clusters in each temper¬ 
ature bin. If the distribution of the luminosity for a given 
temperature is known, we can obtain an average search 
volume by summing Ui„ax(U,L) in a whole range of lu¬ 
minosities weighted with the probability. Assuming that 


luminosities for a given temperature are distributed like 
a Gaussian with a mean given by the power-law function, 
L = AT“, and with a constant standard deviation, in log¬ 
arithmic scale, we can obtain the average search volume 
at a given temperature as 


Unax(T) = 


r w(L,T) 

{log AT°‘ - log Lf 

/ , -exp 

I-oo 



dlogL. (A.6) 


We can then evaluate the temperature function as 

1 


(j){T)AT = 


E 


Knax(Li) ■ 


(A.7) 


Performing Monte-Carlo simulations, we confirmed that 
both estimators, the l/uniax(L,r) method and l/I4iax(T) 
method are consistent and unbiased, that is they repro¬ 
duce the temperature function fed to the simulations in 
the limit of an infinite number of simulated clusters, and 
the mean value of the temperature functions obtained 
from a number of independent measurements with a fi¬ 
nite number of clusters is also equal to the true value in 
the limit of an infinite number of measurements. However, 
the l/Umax(T) method gives smaller variance in a mea¬ 
surement of the temperature function and therefore it is 
more efficient than the l/u„iax(L,T) method. This is well 
demonstrated in Fig. IJ. 

When we apply this analysis method to our cluster 
sample, we further consider the following two points: Our 
first concern is that the flux measurement may have a 
redshift bias such that for clusters with the same true 
luminosity, different luminosities are derived for different 
redshifts. In order to study this effect we produced a sam¬ 
ple of cluster images at various redshifts as observed by 
the LOS'AT'PSPC and determined their fluxes with the 
growth curve analysis method (Bohringer et al. 2000) as 
performed for the actual obse rved d ata of the present sam¬ 
ple by Reiprich & Bohringer ( 2002| ). The simulations show 
that the resulting luminosities are kept constant regardless 
of the redshift — at least in a range that we are consid¬ 
ering here, i.e. z < 0.2, despite the overall signal to noise 
ratio worsening with a dimming of the surface brightness 
that follows (1 -|- z)~'^. 

Secondly, the effect from the measurement errors on 
X-ray fluxes is considered. If there are considerable mea¬ 
surement errors in the flux measurements, a cluster that 
should be in the flux-limited sample could be found to 
be fainter than the flux limit and to be missed from the 
sample, and vice-versa. If the measured flux distribution 
is assumed to be following a Gaussian function with the 
sigma value of df, the probability that a cluster of the lu¬ 
minosity L at the redshift of z is found to be brighter than 
the flux limit fn^ is given as 


W{L,z,h^)= [ 


1 




: exp 


TTCTf 


(/-/o) 


21 


2crf 


df , (A.8) 
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where /o = L/^-kDi^zY. 
modified as 


Therefore, Eq. |A.2| should be 


(T rdV^or 

W(L,T)- 


Setting (Tf to 10 ergs s 


-VE(L,A,/iim)d^ . 


(A.9) 


n“^, a typical flux measure¬ 
ment error around the flux limit, this modihcation brings 
a 2% and 0.5% increase to t^max(-b,T) at 1 and 10 keV, 
respectively, compared with the simple estimate with Eq. 


Although the correction is still insignificant with re¬ 
spect to the Poisson error that each temperature bin would 
have, we use the Eq. A.9 in the derivation of the XTF in 


Sect. 3. Our sample selection criteria specify several pa¬ 
rameters in Eq. A.5, A.6, A.8 and [A.9| as 0=8.139 str, 
/obs.iim = 1-99 X 10“^^ ergs s“^ cm“^ in 0.1-2.4 keV band 
in the observed frame, and then Lobs should be the 0.1- 
2.4 keV band luminosity in the observed frame, which cor¬ 
responds to a cluster rest frame luminosity in the energy 
range from (1 -|- z^ax) 0.1 to (1 -I- Zmax) 2.4 keV. Also, 
from the L — T relation studied in Sect. 2.3 and Sect. 4, 
we estimated the following parameters as log A = 42.15, 
a = 2.47, and criogL = 0.24 for an open universe case, 
and log A = 42.14, a = 2.50, and uiogL = 0.24 for a flat 
universe case. The XTF thus derived is shown in Fig. 


Finally, we should note that Ti in Eq. A.5 is the mea¬ 
sured cluster temperature rather than the “true” temper¬ 
ature. The derived XTF should have been smeared due 
to the temperature measurement error, which brings an¬ 
other systematic error into the resulting XTF. In order to 
investigate the possible systematic errors due to the tem¬ 
perature measurement errors and to confirm our analysis 
method, we performed a Monte-Carlo simulation. With 
an adopted XTF of cx T“®, a number of cluster ob¬ 
servations are simulated taking into account the scatter 
in the L — T relation as well as the flux and tempera¬ 
ture measurement errors, and the XTF w as constructed 
as p erfo rmed with the real data using Eqs. A.5 , A.6, A.8 
and AT. The input XTF and the rederived XTF from the 
simulations show complete agreement except a minor sys¬ 
tematic increase — at most 2% at 10 keV, which is due 
to the temperature measurement error for very hot clus¬ 
ters. These systematics (Eddington bias; Eddington 1944) 
remain in the XTF derived in Sect. 3 (Fig. ^), while, for 
the model htting in Sect. 4, the temperature measurement 
error is taken into account (Eq. 15). 
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Table 1. Summary of the cluster sample 


NAME 

redshift 

nu 

fcTh 

Flux(0.1-2.4keV) 

fix (0.1-2.4keV) {h 

ergs s ^ cm 



(lO^^cm-^) 

(keV) 

(lO'^^ergs s“^ cm“^) 

Open 

Flat 

•PERSEUS 

0.0183 

1.569 


114.33 ± 1.50 

(4.12±0.05)e±44 

(4.19 ± 

0.05)e±44 

•OPHIUCHUS 

0.0280 

2.014 

10.25«j° 

36.25 ± 1.16 

(3.07±0.10)e±44 

(3.14 ± 

0.10)e±44 

COMA 

0.0232 

0.089 

8.07«;f7 

35.13 ±0.79 

(2.04 ± 0.05)e±44 

(2.08 ± 

0.05)e±44 

•A3627 

0.0163 

2.083 


31.60 ±1.57 

(9.04 ± 0.45)e±43 

(9.16 ± 

0.45)e±43 

A3526 

0.0103 

0.825 

Q 5Q+0-05 
0 .uy_Q Q 4 

26.20 ±0.95 

(2.99±0.11)e±43 

(3.02 ± 

0.11)e±43 

•AWM7 

0.0172 

0.921 

Q yf\+0.08 

O. /U_0.05 

16.60 ±0.55 

(5.30±0.18)e±43 

(5.37 ± 

0.18)e±43 

•A2319 

0.0564 

0.877 

8.84t°of4 

12.15 ±0.21 

(4.20 ± 0.07)e±44 

(4.39 ± 

0.07)e±44 

A3571 

0.0397 

0.393 

6.80^°;?! 

12.00 ±0.15 

(2.05 ± 0.03)e±44 

( 2.12 ± 

0.03)e±44 

•TRIANGUL 

0.0510 

1.229 

9.06^;“ 

11.26 ±0.13 

(3.18 ± 0.04)e±44 

(3.31 ± 

0.04)e±44 

A2199 

0.0302 

0.084 

4.28l°;l° 

10.73 ±0.32 

(1.06 ± 0.03)e±44 

(1.09 ± 

0.03)e±44 

•3C129 

0.0223 

6.789 

c ry+O.lG 
-0.15 

10.65 ±0.98 

(5.72 ± 0.53)e±43 

(5.82 ± 

0.54)e±43 

A1060 

0.0114 

0.492 

o 1 r+0.05 
'5.10 _o.o5 

10.04 ± 0.54 

(1.41 ±0.07)e±43 

(1.42 ± 

0.08)e±43 

2A0335 

0.0349 

1.864 

o /- 4 +O.O 9 
'5.0^_0.o8 

8.87 ±0.11 

(1.17±0.01)e±44 

( 1.21 ± 

0.01)e±44 

A0262 

0.0161 

0.552 

9 OK+O-OS 

^.^0_Q QQ 

8.78 ± 0.55 

(2.46±0.15)e±43 

(2.49 ± 

0.16)e±43 

FORNAX 

0.0046 

0.145 

1.00_Q Q7 

8.68 ± 0.66 

(1.97±0.15)e±42 

(1.98 ± 

0.15)e±42 

A0496 

0.0328 

0.568 

4 CQ+O.IO 
^.oy_o.io 

8.20 ± 0.10 

(9.56±0.12)e±43 

(9.81 ± 

0.12)e±43 

A0085 

0.0556 

0.358 

6.5ll°;i® 

7.37 ±0.07 

(2.48 ± 0.03)e±44 

(2.59 ± 

0.03)e±44 

A3667 

0.0560 

0.459 

6 . 281 °;^^ 

7.12 ±0.08 

(2.44 ± 0.03)e±44 

(2.54 ± 

0.03)e±44 

A2029 

0.0767 

0.307 

'7 qq+0-39 
( -yo.Q.se 

6.88 ± 0.07 

(4.43 ± 0.04)e±44 

(4.70 ± 

0.05)e±44 

A3558 

0.0480 

0.363 

c 07 + 0.17 
_o.i5 

6.68 ± 0.05 

(1.68±0.01)e±44 

(1.74 ± 

0.01)e±44 

*S636 

0.0116 

0.642 

2.06^°;°6 

6.65 ± 0.54 

(9.65 ± 0.79)e±42 

(9.75 ± 

0.79)e±42 

A2142 

0.0899 

0.405 

8.46«;g 

6.21 ± 0.09 

(5.52 ± 0.08)e±44 

(5.90 ± 

0.08)e±44 

A1795 

0.0616 

0.120 

^ 17+0.26 
-0.25 

6.21 ±0.03 

(2.58±0.01)e±44 

(2.70 ± 

0.02)e±44 

•PKS0745 

0.1028 

4.349 

6.37«;^i 

6.13 ±0.10 

(7.18±0.11)e±44 

(7.74 ± 

0.12)e±44 

A2256 

0.0601 

0.402 

6.83l°'2? 

6.03 ±0.14 

(2.38 ± 0.05)e±44 

(2.49 ± 

0.06)e±44 

A1367 

0.0216 

0.255 

q p:p:+0.08 

0 . 00 _o.o 8 

5.98 ± 0.07 

(3.02 ± 0.04)e±43 

(3.07 ± 

0.04)e±43 

A3266 

0.0594 

0.148 

7 7 Q+ 0.35 
' • ‘^-0.28 

5.77 ± 0.06 

(2.22 ± 0.02)e±44 

(2.32 ± 

0.03)e±44 

A4038 

0.0283 

0.155 

0 99 + 0.10 
o.^^_ 0 .io 

5.61 ± 0.12 

(4.86±0.10)e±43 

(4.97 ± 

0.10)e±43 

A2147 

0.0351 

0.329 

4 04 + 0.12 
^•9^-0.13 

5.45 ± 0.29 

(7.28 ± 0.39)e±43 

(7.49 ± 

0.40)e±43 

A0401 

0.0748 

1.019 

7 iq+0.28 

‘ -^^-0.24: 

5.26 ± 0.09 

(3.23 ± 0.06)e±44 

(3.41 ± 

0.06)e±44 

N5044 

0.0090 

0.491 

1 99+0.04 

-L-^^_0.04 

5.20 ± 0.04 

(4.54 ± 0.04)e±42 

(4.58 ± 

0.04)e±42 

A0478 

0.0900 

1.527 

fS Q1 +0-40 
'0-y-*--0.36 

5.12 ±0.05 

(4.58 ± 0.04)e±44 

(4.89 ± 

0.05)e±44 

HYDRA-A 

0.0538 

0.486 

q 09 + 0.20 

O.OZ_Q ^7 

4.71 ± 0.04 

(1.49±0.01)e±44 

(1.56 ± 

0.01)e±44 

A2052 

0.0348 

0.290 

0 1 9 + 0.10 
9.1^_0.09 

4.57 ±0.08 

(6.01 ±0.10)e±43 

(6.18 ± 

0.10)e±43 

N1550 

0.0123 

1.159 


4.24 ±0.37 

(6.93±0.61)e±42 

(7.00 ± 

0.62)e±42 

A2063 

0.0354 

0.292 

0 . 00 _Qj ^2 

4.23 ± 0.09 

(5.74±0.12)e±43 

(5.90 ± 

0.12)e±43 

A1644 

0.0474 

0.533 

^4 7n+0-90 
"l- ‘'^-0.70 

4.09 ± 0.34 

(1.00±0.08)e±44 

(1.04 ± 

0.09)e±44 

A0119 

0.0440 

0.310 

<J.Oy_Q 28 

4.05 ± 0.06 

(8.53±0.13)e±43 

(8.82 ± 

0.13)e±43 

•A0644 

0.0704 

0.514 

6-54lg;^^ 

3.97 ±0.07 

(2.16 ± 0.04)e±44 

(2.28 ± 

0.04)e±44 

N4636 

0.0037 

0.175 

0 . 661 H? 

3.97 ±0.47 

(5.84 ± 0.69)e±41 

(5.86 ± 

0.69)e±41 

A3158 

0.0590 

0.106 

r ^1 +0.26 

9'^-*^-0.24 

3.79 ± 0.09 

(1.44±0.04)e±44 

(1.51 ± 

0.04)e±44 

A1736 

0.0461 

0.536 

3.68tgi7 

3.54 ±0.36 

(8.21 ± 0.84)e±43 

(8.51 ± 

0.87)e±43 

A0754 

0.0528 

0.459 

9.00«jt 

3.34 ± 0.09 

(1.01 ±0.03)e±44 

(1.06 ± 

0.03)e±44 

A0399 

0.0715 

1.058 


3.30 ± 0.29 

(1.85±0.16)e±44 

(1.95 ± 

0.17)e±44 

MKW3S 

0.0450 

0.315 

0 4r+0.13 
9.+0_o.io 

3.22 ± 0.05 

(7.13±0.12)e±43 

(7.38 ± 

0.12)e±43 

•A0539 

0.0288 

1.206 

9.U^_0.io 

3.12 ±0.07 

(2.80 ± 0.06)e±43 

(2.87 ± 

0.06)e±43 

EXO0422 

0.0390 

0.640 

tl 9 Qr)+0.90 

^.yU_Q gQ 

3.12 ± 0.32 

(5.16 ± 0.52)e±43 

(5.32 ± 

0.54)e±43 

A4059 

0.0460 

0.110 

0 Q4+0.15 
9.y^-o.i5 

3.10 ±0.07 

(7.16±0.16)e±43 

(7.42 ± 

0.16)e±43 

A3581 

0.0214 

0.426 

1 Qq+0.04 

T.OO_q Q 2 

3.08 ± 0.16 

(1.53±0.08)e±43 

(1.55 ± 

0.08)e±43 

A3112 

0.0750 

0.253 

'^-0.25 

3.07 ±0.06 

(1.90±0.03)e±44 

( 2.01 ± 

0.04)e±44 

A0576 

0.0381 

0.569 

q oq+ 0 .i 6 
O.00_o.i5 

2.99 ± 0.33 

(4.73 ± 0.53)e±43 

(4.88 ± 

0.54)e±43 

A3562 

0.0499 

0.391 


2.89 ± 0.04 

(7.87±0.12)e±43 

(8.18 ± 

0.13)e±43 

A2204 

0.1523 

0.594 

6.38^°;^^ 

2.73 ± 0.07 

(7.13±0.19)e±44 

(7.92 ± 

0.21)e±44 

A0400 

0.0240 

0.938 

2.43^°;!^ 

2.64 ± 0.05 

(1.65±0.03)e±43 

( 1.68 ± 

0.03)e±43 
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Table 1. continued 


NAME 

redshift 

riH 

fcTh 

Flux(0.1-2.4keV) 

“■Lx (0.1-2.4keV) {h 

^ ergs s ^ cm 



( 10 ^^cm~ 2 ) 

(keV) 

(lO'^^ergs s“^ cm~^) 

Open 

Flat 

A2065 

0.0721 

0.284 

c q7-l-u.d4 
O.Oi _0.30 

2.55 ±0.26 

(1.46±0.15)e±44 

(1.54 ± 

0.15)e±44 

A1651 

0.0860 

0.171 

6 

'->•^^-0.41 

2.52 ± 0.05 

(2.06 ± 0.04)e±44 

(2.19 ± 

0.04)e±44 

A2589 

0.0416 

0.439 

q qq+0.13 

o.oo_Q 

2.52 ± 0.05 

(4.75 ± 0.10)e±43 

(4.91 ± 

0.11)e±43 

MKW 8 

0.0270 

0.260 

q 9Q+0.23 

O-^y_ 0.22 

2.51 ±0.35 

(1.99±0.27)e±43 

(2.03 ± 

0.28)e±43 

A2657 

0.0404 

0.527 

q t::q+ 0.12 

o.oo_Q 9^2 

2.49 ± 0.04 

(4.44 ± 0.06)e±43 

(4.58 ± 

0.07)e±43 

A3376 

0.0455 

0.501 

4.43l°j9 

2.44 ± 0.06 

(5.52 ± 0.13)e±43 

(5.72 ± 

0.14)e±43 

SllOl 

0.0580 

0.185 

^2.60lg;®° 

2.44 ± 0.04 

(9.05 ± 0.14)e±43 

(9.46 ± 

0.14)e±43 

A1650 

0.0845 

0.154 

5.68l°;l? 

2.43 ± 0.26 

(1.92 ± 0.21)e±44 

(2.04 ± 

0.22)e±44 

A2634 

0.0312 

0.517 

3-45l°;l® 

2.38 ± 0.06 

(2.51 ± 0.06)e±43 

(2.58 ± 

0.07)e±43 

A3391 

0.0531 

0.542 

5.89l°;« 

2.20 ± 0.07 

(6.78 ± 0.21)e±43 

(7.06 ± 

0.22)e±43 

A2597 

0.0852 

0.250 

^•^'^-0.41 

2.20 ± 0.04 

(1.77 ± 0.04)e±44 

( 1.88 ± 

0.04)e±44 

ZwC11215 

0.0750 

0.164 

^6.30j;i-^^ 

2.17 ±0.05 

(1.34 ± 0.03)e±44 

(1.42 ± 

0.03)e±44 

A2244 

0.0970 

0.207 

r yy+0.61 

• • -0.44 

2.10 ±0.07 

(2.19 ± 0.08)e±44 

(2.35 ± 

0.08)e±44 

A0133 

0.0569 

0.160 

q Q 7 + 0.28 
-0.27 

2.06 ± 0.03 

(7.33 ± 0.10)e±43 

(7.66 ± 

0.11)e±43 

A2163 

0.2010 

1.227 

lU.JO_o .68 

2.04 ± 0.05 

(9.28 ± 0.23)e±44 

(10.59 ± 0.27)e±44 

A2255 

0.0800 

0.251 

r 09+0.40 

O.y^_0.26 

2.02 ± 0.04 

(1.42 ± 0.03)e±44 

(1.51 ± 

0.03)e±44 

IIIZw54 

0.0311 

1.668 


2.01 ± 0.25 

(2.11 ± 0.27)e±43 

(2.16 ± 

0.27)e±43 

A3395S 

0.0498 

0.849 

c p:p:+0.89 
J.JJ_0.65 

2.01 ± 0.13 

(5.42 ± 0.34)e±43 

(5.63 ± 

0.35)e±43 

N507 

0.0165 

0.525 

1 ^q+0.04 
J-."±U_o.o7 

2.00 ± 0.04 

(5.88 ± 0.12)e±42 

(5.96 ± 

0.12)e±42 

MKW4 

0.0200 

0.186 

-l--o4_0 03 

2.00 ± 0.05 

(8.65 ± 0.24)e±42 

(8.79 ± 

0.24)e±42 

UGC03957 

0.0340 

0.459 

b q 1 0+1.36 

J-iy_0.95 

Uq 91 +1.48 

O-^l_ 1.02 

1>7 t::t::+3.21 
/ •JJ_2.25 

1.98 ±0.19 

(2.49 ± 0.24)e±43 

(2.56 ± 

0.25)e±43 

A3827 

0.0980 

0.284 

1.98 ±0.19 

(2.10 ± 0.20)e±44 

(2.25 ± 

0.22)e±44 

“7.66l|®^ 

A3822 

0.0760 

0.212 

c 19+0.43 
J-l^_0.31 

1.98 ± 0.24 

(1.26±0.15)e±44 

(1.33 ± 

0.16)e±44 

IIZwl08 

0.0494 

0.663 

^4.25l}'®i 

»4.28ll;®« 

1.90 ±0.23 

(5.06 ± 0.61)e±43 

(5.26 ± 

0.63)e±43 

M49 

0.0044 

0.159 

1 qq + 0.03 
i.oo_0.03 

1.89 ± 0.03 

(3.93 ± 0.07)e±41 

(3.94 ± 

0.07)e±41 

ZwC11742 

0.0757 

0.356 

“6.05111? 

1.87 ±0.13 

(1.18 ± 0.08)e±44 

(1.25 ± 

0.09)e±44 

S405 

0.0613 

0.765 

'’^■98tlil 

“5.02l?j2 

1.82 ± 0.24 

(7.48± 1.00)e±43 

(7.84 ± 

1.05)e±43 

A3532 

0.0539 

0.596 

4-411^18 

1.78 ± 0.05 

(5.65 ± 0.17)e±43 

(5.88 ± 

0.17)e±43 

A3695 

0.0890 

0.356 

“6.67l?J^ 

“6.76+?+ 

1.76 ± 0.26 

(1.54±0.23)e±44 

(1.65 ± 

0.25)e±44 

HCG94 

0.0417 

0.455 

3.30+;l^ 

1.72 ± 0.03 

(3.25 ± 0.05)e±43 

(3.36 ± 

0.06)e±43 

A3528S 

0.0551 

0.610 

4.60++ 

1.70 ± 0.04 

(5.65 ± 0.15)e±43 

(5.89 ± 

0.15)e±43 

S540 

0.0358 

0.353 

b q 07+1-30 
o.U< _o.9i 

“3.09++ 

1.62 ± 0.13 

(2.26 ± 0.18)e±43 

(2.32 ± 

0.19)e±43 

A2877 

0.0241 

0.210 

“3.50+??? 

1.61 ±0.03 

(1.01±0.02)e±43 

(1.03 ± 

0.02)e±43 

A3395n 

0.0498 

0.542 

cr -,-,+0.47 
^•-1-1-0.43 

1.54 ± 0.10 

(4.16 ± 0.27)e±43 

(4.33 ± 

0.28)e±43 

A2151w 

0.0369 

0.336 

Z.JO_o.20 

1.53 ± 0.05 

(2.26 ± 0.07)e±43 

(2.33 ± 

0.07)e±43 

A3560 

0.0495 

0.392 

“3.871?-?? 

“3.901?-?? 

1.50 ± 0.06 

(4.02 ± 0.17)e±43 

(4.17 ± 

0.17)e±43 

A2734 

0.0620 

0.184 

5.07+;?? 

1.47 ±0.06 

(6.19 ± 0.25)e±43 

(6.49 ± 

0.26)e±43 

A0548e 

0.0410 

0.188 

9 00+0.17 
^.yo_o.i5 

1.46 ± 0.04 

(2.68 ± 0.08)e±43 

(2.77 ± 

0.08)e±43 

A1689 

0.1840 

0.180 

8.581?;?? 

1.45 ± 0.03 

(5.53 ± 0.10)e±44 

(6.25 ± 

0.11)e±44 

A1914 

0.1712 

0.097 

8.41+;?? 

1.45 ± 0.03 

(4.77±0.11)e±44 

(5.36 ± 

0.12)e±44 

RXJ2344 

0.0786 

0.354 

“5.451?-?? 

Ur r9+2.55 
J*J-^_1 74 

1.37 ±0.03 

(9.34±0.21)e±43 

(9.91 ± 

0.23)e±43 

A3921 

0.0936 

0.280 

5.39+;?? 

1.31 ± 0.04 

(1.27±0.03)e±44 

(1.36 ± 

0.04)e±44 

A1413 

0.1427 

0.162 

6.56+;?? 

1.28 ± 0.03 

(2.92 ± 0.08)e±44 

(3.23 ± 

0.09)e±44 
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Table 1. continued 


NAME 

redshift 

riH 

fcTh 

Flux(0.1-2.4keV) 

1'Lx (0.1-2.4keV) {h 

^ ergs s ° cm ‘‘) 



(lO^^cm-2) 

(keV) 

(10~“ergs s"^ cm"^) 

Open 

Flat 

N5813 

0.0064 

0.419 

0.761°:;“ 

1.27 ±0.14 

(5.62 ± 0.62)e±41 

(5.66±0.62)e±41 

A1775 

0.0757 

0.100 

3.66l°:l^ 

1.26 ±0.04 

(8.02 ± 0.23)e±43 

(8.49 ± 0.25)e±43 

A3528n 

0.0540 

0.610 

4 7Q+0-50 
' ^-0.44 

1.26 ±0.05 

(4.00 ± 0.15)e±43 

(4.17±0.16)e±43 

A1800 

0.0748 

0.118 

'4.95i;:“ 

1.20 ±0.15 

(7.39 ± 0.95)e±43 

(7.82 ± 1.01)e±43 

A3888 

0.1510 

0.120 

^8.46l2'53 

»8.68l^:?i 

1.09 ±0.04 

(2.77±0.11)e±44 

(3.08±0.12)e±44 

A3530 

0.0544 

0.600 

4t.UO_o.30 

0.96 ±0.04 

(3.11±0.15)e±43 

(3.25±0.15)e±43 

N5846 

0.0061 

0.425 

U.D4_0.03 

0.83 ±0.03 

(3.35 ± 0.12)e±41 

(3.36±0.12)e±41 

N499 

0.0147 

0.525 

U.DU_0.03 

0.45 ±0.02 

(1.07±0.04)e±42 

(1.08±0.05)e±42 

A0548w 

0.0424 

0.179 

-L.U^_0.49 

-L.UO o.53 

0.25 ±0.02 

(4.84 ± 0.43)e±42 

(5.03 ± 0.45)e±42 


: The Galactic latitude is lower than 20°. 

d : The values are calculated for an open (flm,o = 0.2, flA.o = 0) and a flat (flm.o = 0.2, flA.o = 0.8) universe. 
: The temperature is estimated with the L — T relation assuming the open universe. 

: The temperature is estimated with the L — T relation assuming the flat universe. 

H) 


: The temperature is derived with non-ASCA spectroscopy. Einstein MFC for A2877 and A1644 (David et al. 1993), EXOSAT 
LE-I-ME for EXO0422 (Edge & Stewart 1991), and XMM-Newton EPIC for SllOl (Kaastra et al. 2001). 


Table 3. Fitting Results 


Model XTF 

^m ,0 

(78 

A 

a 

L7log L 

PS(Open) 

PS(Open,r > 3keV) 

FE(Open,s=0) 

FE(Open,s=l) 

PS (Flat) 

PS(Flat,r > 3keV) 

FE(Flat,s=0) 

FE(Flat,s=l) 

n 1 q'FU'jo 

U.io_o.o5 

U.^<-)_0.09 
f) 1 c:+0.07 

U.iO_o.o5 

^•-*--*--0.05 

0 iq +0 08 

U.iy_0.05 

0-26l°:J= 

o.i8l°:°5 

0.141^5 

0.96:°:“ 

0.891°:;? 

0.98in? 

0.681^3 

i.02i°:;? 

0.941°:;^ 

i.i4i°:;? 

0.89in^ 

42 1 

4tZ.iO_o.20 

4±z.iy_o 28 

42 14+0-1'^ 
4±Z.14i_o 20 

4tZ.lZ_o 20 

42 1 c^+017 

4tZ.iO_o.20 

42 IQ+O-27 
4iZ.iy_o 29 

42.i5i°:;j 

42.141°:;; 

2 47 +u.^t. 

-0.24 

9 A 9+0-38 

^•4tZ_o 37 

2 4Q+°-20 

^•4ty_o 25 

9 +0-27 

^•’->i_o.25 

2 4Q+°-20 
^•4ty_o 25 

2 44+°-3^ 
^•^^-0.38 

2 4Q+°-20 

^•4ty_o 25 

9 +0-26 
^•’->i o.25 

U.Z4i_o 04 
r\ qq + 0.05 
u.zz_o 03 
n 24+0 04 

U.Z4i_o 04 
n 24+0-04 
U.Z4i_o 04 

n 24+0-04 
U.Z4i_o 03 

90+0.04 

U.z,O_0.04 

0 24+0-04 
U.Z4i_o 03 

0 24+0-04 
U.Z4i o 03 


The error ranges corresponding to 90% confidence for one parameter of interest. 













